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ABSTRACT 

I model profiles of the [Ne n] forbidden emission line at 12.8 lyum, emitted by photoevaporative 
winds from discs around young, solar-mass stars. The predicted line luminosities (~ lO^^L©) 
are consistent with recent data, and the line profiles vary significantly with disc inclination. 
Edge-on discs show broad (30^0kms ') double-peaked profiles, due to the rotation of the 
disc, while in face-on discs the structure of the wind results in a narrower line 10km s ') 
and a significant blue-shift (5-lOkms"'). These results suggest that observations of [Nen] 
line profiles can provide a direct test of models of protoplanetary disc photoevaporation. 
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1 INTRODUCTION 

The manner in which gas is removed from discs around young 
stars has critical consequences for theories of planet formation, as 
, the formation of gas-giant planets must precede the dissipation of 
' gas discs. Currently, models of protoplanetary disc evolution sug- 
I gest that the most of the gas is removed by a combination of vis- 

• cous accretion through the disc and photoevaporation by irr a diatio n 
' from the central sta r [e.g.. Clarke. Gendrin & Sotomavoil 1 I2OOII) : 
, [Alexander. Clarke~& Pringle (2006bi): see a l so the recent reviews 
■ bv iDullemond et al.l ( |2007|) and lAlexandeil ( |2008|) 1. Such models 

agree reasonably well with current data, but the ionizing luminosi- 
ties of T Tauri stars (re quired to drive the photoevaporation ) are 
very poorly constrained. Alexander. Clarke & Pringlg ( l2005h used 
archival data to estimate the chromospheric ionizing fluxes from a 
, small sample of bright T Tauri stars, and found values in the range 

• ~ lO'^^-lO'*^^ ionizing photons per second. However, recent data 
suggests that this may ov er-estimate mo re typical values of O by 

1 1-2 orders of magnitude jHercze Photoevaporative flows 

have been obse rved clearly in ul t racom p act Hii regions around 
massive stars ( iHollenbach et alj 1 1994 iLugo. Lizano. & Garavl 

l2004h . and in the externally illuminated "proplyds" in the Orion 
nebula (e.g. |jphnstone, Hollenbach & Ballv 1998). However, the 
evidence for or against the existence of central star-driven pho- 
toevapor ative wind s from discs around solar-mass stars is rather 
tenuous. iFont et al] ( |2004|) modelled the profiles of optical forbid- 
den emission lines and found reasonable agreement with both the 
line fluxes and profiles obtained in t he spectroscopic survey of 

iHartigan. Edwards & GhandouJ Il995h . However, some discrepan- 
cies were also present (notably between the predicted and observed 
fluxes from [Oi]), and no further studies of such emission lines ex- 
ist. 

For several reasons, the [Nen] fine-structure line at 12.81;/m 
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offers perhaps the best opportunity to confirm or deny the existence 
of a slow (~ 10 km s"') ionized wind from a T Tauri disc. The high 
ionization potential of neon (21.56eV) means that Ne^ only exists 
close to T Tauri stars in photoionized gas, so forbidden emission 
lines from neon ions are unlikely to arise elsewhere in the T Tauri 
system. The 12.81yum line also falls in an atmospheric window (un- 
like the [Nein] line at 15.55;um), and can therefore be observed 
from the ground at echelle resolution. Moreover, the continuum 
emission from the star-disc system is orders of magnitude weaker 
in the mid-infrared than in the optical, so the line-to-continuum ra- 
tio for [Ne n] can often be higher than for optical forbidden lines 
such as [Nil] or [Sii]. In addition, as will be shown, the critical 
density of the [Nen] line is such that most of the emission arises in 
the "launching region" of the photoevaporative wind, making it an 
excellent tracer of disc photoevaporation. 

Recently the Spitzer Space Telescope has detected the [Ne n] 
12.81//m line towards more than 20 y o ung, solar-mass star s 
jPascucci et al] |2007| : iLahuis et all |2007| : lEspaillat et all |2007|) . 
and in most cases it is thought to arise in the star-disc system. 
iHerczeg et alj ( |2007|) observed the [Nen] line from the nearby 
source TW Hya at high resolution (A/AA » 30,000), and found 
that the line (with width 2 1 ± 4 km s" ' ) was narrower than lines that 
originate in the accretion funnel, but significantly broader than lines 
emitted from the disc at larger radii: this strongly suggests that the 
line or igina tes in a hot disc atm o sphere . iGlassgold. Naiita & Igeal 
1 I2OO7I) and iGorti & HoUenbacH ( l2008h showed that X-ray heat- 
ing of the disc can reproduce the observed line fluxes, but a pho- 
toevaporative disc wind should give rise to a similar line flux 
jPascucci et ai]|2007l : iGorti & Hollenbacril200^ . While both ion- 
ization mechanisms (X-ray and UV) can likely explain the observed 
line fluxes, the velocity structure of a photoevaporative wind should 
result in a line profile that is distinct from that produced by X-ray 
irradiation of a bound disc atmosphere. In this Letter I model the 
profiles of the [Nen] 12.8 1/jm line emitted by photoevaporative 
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winds, and find that high-resolution spectroscopy can provide an 
unambiguous means of detecting such winds in T Tauri systems. 



2 MODELS 

2.1 Hydrodynamic models 

In order to construct velocity profiles it is first necessary to 
compute the hydrodynamic structure of the wind. Here two dif- 
ferent cases are considered: the standard photoevaporative wind 
jHoUenbach et alj|l994l : [pont et alj|2004h and the "direct" wind, 
applicable to a disc with an inner hole that is optically thin to ioniz- 
ing photons (Alexander, Clark e^ Pring le 2006aJ. I use the zeus2d 



hydrodynamics code Jstone & NormanI 



199: 



to model the winds, 

followi ng the approaches of lFont et al.l ( |2004) and lAlexander et al] 
l l2006al) respectively. I assume azimuthal and midplane symmetry, 
and therefore model the disc using a polar [(r, 6*)] grid covering the 
range 9 = [0, ;r/2]3. The grid cells are logarithmically-spaced in r 
and linearly-spaced in 6, with the numbers of radial (A', ) and po- 
lar (Ng) cells chosen so that the grid cells a re approximately square 
throughout (i.e. Ar = rAS, see e.g., lBate"e t al. 2002). The rotation 
option in zeus2d, which introduces a centrifugal pseudo-force, is 
turned on, and accelerations due to gravity are evaluated using only 
a point mass at the origin. I adopt the van Leer (second-order) in- 
terpolation scheme, and the standard von Neumann & Richtmyer 
form for the artificial viscosity (with ^visc = 2.0). 

2.1.1 Standard case 

In the standard case th e wind models fol low the "photoevaporative 
disc wind" mod els oflFontet"^ i2004[)- These models make use 



19941 specifically the "weak- 



of the results of iHoUenbach et ^ 
wind" case), who used detailed radiative transfer calculations to 
determine the structure of the ionized disc atmosphere. The atmo- 
sphere is isothermal at lO^'K, with a sound speed c, = lOkms"'. 
The critical length-scale R^, known as the gravitational radius, is 
found by equating the sound speed of the ionized gas with the local 
(Keplerian) orbital speed, so 
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M. 



AU, 
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where M, is the mass of the central star The (number) density at 
the base of the ionized atmosphere, no(R), determines the structure 
of the wind, and the value at R„ is 
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atomic hydrogen at 10''K ('Cox"2000^, and the constant C was de- 
termined by the numerical calculations of Hollenb ach et alj ( ll994h . 
Inside /? » the base density scales as iJ"^'^, while outside i?g it scales 
as J?~^^^. lFont et al.l j2004l) adopt the following fitting form for the 
base density, which varies smoothly between the two power-laws: 



no(R) 



[(R/R„y^''- + {R/Rgy^'^ 



1/5 



(3) 



The disc atmosphere is given an isothermal equation of state, and 
the outer and inner radial boundaries are set as outflow boundaries. 



Lower-case r denotes spherical radius; upper-case R cylindrical radius. 




Figure 1. Density and velocity structure of the steady-state photoevapora- 
tive flow, z = corresponds to the base of the disc atmosphere; the under- 
lying disc is not shown. The plot shows the "launching region" of the wind; 
the grid extends to r = 20Rg . The density is shown on the colour-scale, with 
velocity vectors plotted at regular intervals. The white line shows the sonic 
surface, where the velocity in the (r, 9) plane is equal to the sound speed. 
The black contour indicates where the local density is equal to the critical 
density for the [Neii] line, for parameters <D = 10'" s"' and M, = IMq. 



The inner polar boundary (the z-axis) is a reflecting boundary, but 
the flow is not at all sensitive to the inner boundary conditions. 
The base of the flow (i.e. the disc midplane) is held constant at 
every time-step, using the density profile no(R) from Equation |2] 
[The base of the flow is actually found several disc scale-heights 
above the midplane (Hollenbach et al. 1994), but additional cal- 
culations show that this approximation has a negligible effect on 
the line profiles.] The flow is thus assumed to be recombination- 
dominated (with negligible advection across the ionization front): 
previous studies suggest that this approximation is valid (e.g. 
[Alexander e t al. 2006a). The base cells are given a Keplerian veloc- 
ity in the orbital direction, and zero velocity in the radial direction. 
The polar (vertical) velocity out of the base cells is not prescribed, 
but instead computed hydrodynamical ly: typically the launch ve- 
locity is 0.3-0.4C, dPont et al.P 2004). The models are computed 
in dimensionless units: the unit of length is R„, the unit of time is 
the orbital period at R„, and the density is normalized so that = 1. 
Initially the grid is filled with a uniform density of lO^^ng, and the 
model is evolved forward in time until a steady-state is reached. 
In order to compute line profiles, the resulting density and velocity 
structures are scaled to physical units by choosing the parameters 
O and M,. 

I use a grid with r,„ = 0.03Rg, r„u, = 20Rg, A', = 246 and 
Ng = 59. As one would expect there is excellent agreement with the 
results of iFont et al.l ( I2004D . with the location of the sonic surface 
and the integrated mass-loss rate agreeing with their result to better 
than 1% accuracy. Choosing different locations for the boundaries 
(rin = O.OlRg; r^ut = 50Rg) has a negligible effect on the structure 
of the flow and affects the computed line fluxes and profiles only 
at the percent level. A numerical convergence test (using twice the 
grid resolution in both dimensions) suggests that the flow solution 
is correct to 1-2% accuracy. A steady-state is reached after t = 
5-10 time units; the subsequent line profile modelling makes use 
of the density and velocity fields at ? = 40. The steady-state flow 
solution is shown in Fig[T] 
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2.1.2 Discs with holes 

For the case of a disc with an optically thin inner cavity, the model 
of [Alexander et al. (2006a) is used. In this case the calculation is 
not scale-free (because of the ionization/recombination balance cal- 
culation), so physical units are used. In addition, because the pres- 
sure scale-height of the cold disc must be resolved into several grid 
cells, somewhat higher resolution is required than in the standard 
case (where the underlying cold disc is treated as a boundary condi- 
tion). I adopt model parameters of M, = IMq and (t> = lO'*' s"' , and 
grid parameters r,„ = 5AU, r^^t = 200AU, N,. = 360 and Ng = 153. 
The outer radius is chosen to be approximately equal to that used 
in the standard case when scaled to physical units (for M, = IMq), 
but the results are not sensitive to the exact location of the bound- 
ary. The initial disc follows an power-law surface density profile 
Z oc R ', truncated exponentially at some inner radius, and is lo- 
cally isothermal in the vertical direction with aspect H/R = 0.05. 
The total disc mass is O.OOIMq (consistent with the low disc mass 
expected during photoevaporative clearing), but the resulting line 
profiles are not sensitive to this mass. The model is allowed to 
evolve to a "quasi-steady" stat^ before computing line profiles, 
and considers a inner hole of radius 9AU (-Rg). The quasi-steady 
state is reached after 1-2 outer orbital periods; the line profiles are 
calculated from the density and velocity fields at ? = 4000yr. 



2.2 Line profiles 

In order to compute line profiles from the hydrodynamic models, it 
is first necessary to construct three-dimensional density and veloc- 
ity fields from the two-dimensional simulations. I assume reflective 
symmetry at the disc midplane in order to extend the coordinate 
range lo 6 = [0,n], and azimuthal symmetry around the polar axis. 
The resulting grid has N,- x 2Ng x cells: an azimuthal resolution 
of A'^ = 120 is adopted, but the exact value of has no significant 
effect on the results. 

The line luminosity L at a given velocity v is computed as 



L(V): 



I 



V27r' 



TTVth 



|exp(-I^ 



A^Ne Xll 



xn,(r)P„A„i hv.idV, 



(4) 



with the integral evaluated by direct summation over the entire ex- 
tent of the grid. Here vios(r) is the line-of-sight component of the 
local gas velocity vector, computed for a specified inclination an- 
gle /. (( = ;r/2 corresponds to a disc viewed edge-on; i = face-on.) 
ng(r) is the local number density of the gas (hydrogen is assumed 
to be fully ionized), and hv„i is the energy of the emitted photons. 
The Doppler broadening term depends on the thermal velocity of 
the emitting atoms 



V,h = Cs V^Ww^Ne ' 



(5) 



where and ntf^^ are the masses of hydrogen and neon atoms re- 
spectively. I adopt the standard solar value for the abundance of 
neon, A/j^e = 10 x lO"'', a nd the Einstein coefficient of the transi - 
tionisA,,, = 8.39 x 10"^ s"' ( lMendozalll983l ; lGlassgold et al.l2007h . 



^ The flow in this case is only "quasi-steady" because the mass reservoir is 
finite. This state is reached in < lO'yr, which is much less than the time- 
scale for disc clearing (~ lO^yr). As a result, aside from initial transients 
the flow solution does not change significantly over the duration of the sim- 
ulations. 



Xll is the fraction of neon that exists as Ne^, discussed below. Fol- 
lowing Glassgold et al. (2007i), the excitation fraction of the upper 
state P„ is computed as 

1 



P. 



(6) 



2C„,exp(-r„,/r)-H I 

where the excitation temperature r„; = 1122.8K, and the gas tem- 
perature T = 10, OOOK. C„; expresses the departure from a thermal 
level population, and is defined as 



C,„ = 1 + 



riAr) 



(7) 



where the critical density n^T = 5 x lO^cm"^. I assume that the 
wind is optically thin to the emitted line, but that the disc absorbs 
100% of the line flux if the midplane intercepts the line-of-sight to 
the observer. This approach neglects line emission from partially 
ionized gas in the ionization front, but analytic estimates suggest 
this contributes to the total line luminosity only at the percent level. 

The factor Xn is determined by ionization/recombination bal- 
ance. The first ionization potential of atomic neon is 21.56eV, so in 
the lO^'K disc atmosphere collisional ionization of neon is negligi- 
ble. Instead, neon is photo-ionized, either by UV photons with ener- 
gies greater than 21.56eV or by high-energy X-rays (via the Auger 
mechanism). The low line-of-sight column density in the wind re- 
gion (10'''-10'*'cm"^) suggests that UV ionization is more efficient 
than ionization by X-rays, and I neglect X-ray ionization here. The 
second ionization potential of neon is 41.0eV, so the relative abun- 
dances of Ne, Ne+ and Ne^^ are determined by the spectral slope of 
the incident UV radiation field. Little is known about the details of 
the incident radiation field, but comparison to Galactic H ii regions 
suggests that ^"11 likely lies in the range 0.1-1.0 (e.g. iRubin et al.l 
1I99 li) . In the models almost all the ionized gas is below the critical 
density, so the line flux per unit volume is simply proportional to 
Xll. Moreover, the bulk of the line flux arises from a fairly limited 
range in radius (^ 0. l-2.0^g; see discussion in Section |3] below), 
and Xll is unlikely to vary significantly over this range. Ionization 
balance therefore has a negligible effect on the shape of the line 
profile, but the integrated line flux scales linearly with Xu. In the 
absence of any knowledge of the incident spectrum of radiation I 
adopt a constant value of Xu = 1/3, noting that this introduces an 
uncertainly of a factor of ~ 3 in the derived line luminosities. 



3 RESULTS 

The results of the line-profile modelling are shown in Table[T] Two 
disc models are presented (standard, 9AU hole) for model param- 
eters M, = IMq and an O = lO^'s"', and for each model the line 
profile was computed for inclinations of i = (face-on), i = n/4 
and I = 7r/2 (edge-on). In addition, "realistic" line profiles, at 
the typical resolution of an echelle spectrograph, were computed 
by convolving the line profiles with a Gaussian profile with half- 
width cr = 5kms"': the effective resolution of these profiles is 
A/AA =30,000. 

The line profiles for the standard case are shown in Fig|2l 
for clarity, only the edge-on and face-on cases are shown. When 
viewed edge-on the line is broad (30-40 km s"') and double- 
peaked, while face-on inclination results in a line that is narrower 
(lOkms"') and slightly blue-shifted (6-7 km s"'). At resolution 
A/AA = 30,000 the blue-shift is comparable for all inclinations 
i < n/4 (see Table [T]l, as any double-peaked structure in the line 
is not well-resolved. These line profiles can be easily understood 
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Inclination 


Line flux 


Peak FWHM 


Peak 


FWHM 




10-'' Lo 


Model 


Xjl^X 


= 30,000 






Standard model 






n/2 


3.0 


±10.6 33.5 


±7.6 


38.1 




1.5 


-12.1/+0.9 26.3 


-7.4 


28.7 





1.5 


-6.4 9.5 


-6.6 


16.8 




Disc with R„-sized inner hole 






nil 


5.4 


±3.1 17.1 


0.0 


19.3 


jr/4 


2.7 


-4.8 13.3 


-6.3 


18.1 





2.7 


-7.9 10.5 


-8.9 


17.0 



Table 1. Details of the line profiles shown in Figs|2]&[3] computed for a 
stellar mass M» = IM© and an ionizing flux <S) = 10^' s"'. The columns 
labelled "Peak" show the velocities at which the line fluxes are maximum, 
while the full widths of the hnes at half the maximum flux are labelled 
"FWHM". "Model" denotes the values for the theoretical line profiles, 
while "A/AA = 30,000" indicates the values when observed at the typical 
resolution of an echelle spectrograph. Velocities with positive and nega- 
tive values indicate that the profile is double-peaked. All velocities are in 
km s"' ; the ambient disc gas has zero velocity. 
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Figure 3. As Fig|2] but for a disc with a 9AU inner hole. The face-on line 
profile is similar to that seen in the standard case, but in the edge-on case the 
absence of the inner disc results in a narrower line, with a less pronounced 
double-peaked profile. 
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Figure 2. Line profiles for the photoevaporative disc wind. The black lines 
show the profiles when the disc is viewed edge-on (r = ;r/2), the red lines 
the profiles for a face-on disc ((' = 0). The dashed lines show the profiles 
degraded to typical echelle resolution (/i/A/1 = 30,000). The flux per unit 
velocity interval is normalised to the peak value in the edge-on case. When 
viewed edge-on the fine is broad and double-peaked, reflecting the large, 
10-20 km s"', orbital velocity in the emitting region. When viewed face-on 
the disc rotation is perpendicular to the line-of-sight, and instead the vertical 
component of the wind results in the line being blue-shifted by 5-10 km s"' . 



by considering the flow structure. At large radii the wind is essen- 
tially isothermal and spherically symmetric, and is analogous to a 
Parker wind solution. The streamlines are close to radial, and the 
density along streamlines drops as oc r"^. As the density is sub- 
critical the line flux per unit volume scales oc oc r"'', so the total 
line luminosity (integrated over volume) is dominated by the emis- 
sion fro m small radii. In spection of the mass-loss profile (see, e.g., 
Fig.7 of lFont et al.l2004l) shows that > 80% of the mass-loss comes 
from < IRg, and the density structure is such that the critical den- 
sity is reached only at radii < O.li^u (see Fig[T). This highlights 
why the [Neii] 12.8 lyum line is an ideal probe of disc photoevapo- 
ration: lines with higher critical densities will likely be dominated 



by emission from high density regions very close to the star, while 
lines with lower critical densities are insensitive to the launching 
region of the wind. In the launching region the rotation speed of 
the disc is greater than the wind or thermal speeds (as vk = t's at 
R = R,,), so when viewed edge-on the line profile is dominated by 
Keplerian rotation and shows a pronounced double-peaked struc- 
ture. When viewed face-on the rotation is entirely in the plane of 
the sky, and instead the vertical component of the wind velocity 
results in a blue-shifted profile. In this case the profile is slightly 
asymmetric, with a blue "tail" due to high-velocity, low-density 
gas at large r. In all cases the line-widths are significantly larger 
than expected from thermal broadening alone. This occurs because 
the streamlines are almost radial at large radius. Consequently the 
dispersion in the line-of-sight component of the flow velocity is al- 
ways larger than the thermal line-width, resulting in broader than 
thermal line- widths even for face-on discs. 

In the case of a disc with a cleared central hole, the line profile 
when viewed edge-on is markedly different from the standard case 
(see Fig[3}. The removal of gas with high rotation speeds results in 
a profile that is much less strongly double-peaked, and significantly 
narrower (15-20 kms"'). Still larger hole sizes result in even less 
pronounced double-peaked profiles, and for hole sizes > 5Rg there 
is no evidence of a double-peaked profile. The profile in the face-on 
case is similar to that seen in the standard model: fairly narrow (^ 
lOkms"'), slightly asymmetric, and blue- shifted (7~9kms"'), and 
this general shape is seen irrespective of the size of the central hole. 
The integrated line flux is somewhat larger than in the standard case 
(by a factor of =^ 2), due to the higher density in the wind launching 
region. 

The predicted line luminosities are a few xlO^^L©. The ex- 
act line fluxes are likely only accurate to within a factor of ~ 3, 
as discussed in Section [2!2l above, but the predicted line strengths 
are comparable to t hose predicted by previous calculations, both 
from UV ionization ( Gorti & Hollenbactill2008l) and X-ray ioniza- 
tion jGlassgold et"ai]|2007t) . The emitted line fluxes scale approxi- 
mately linearly with <I>, and the shape of the line profile is largely 
independent of the value of O. This can be understood from Equa- 
tions O&IH the line flux scales oc as long as <c and the 
density of the ionized gas scales oc O''^. We see from Fig[T]that 
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almost all of the gas in the emitting region is below the critical den- 
sity for [Ne ii], and additional calculations suggest that this holds as 
long as O < IC'^s"'. Higher still ionizing fluxes will result in the 
critical density being exceeded, making the [Nen] line less sensi- 
tive to the launching region of the wind, but this is unlikely to be 
relevant to discs around T Tauri stars. 



4 DISCUSSION 

In general, t hese results com pare favourably with the previous 
modelling o f iFont et alj ( [2004 . These authors modelled the emis- 
sion from optical forbidden lines ([Sii] and [Nii]), and for face- 
on inclinations found that the line profiles were blue-shifted by 
^ lOkms"'. The blue-shift decreased with increasing inclination 
angle, and in the edge-on discs the lines were centred on zero ve- 
locity. They did not, however, find any double-peaked lines, even in 
the edge-on case. The reasons for this are not entirely clear, but it is 
likely attributable to the different critical densities of the lines con- 
sidered. The critical densities of the [S ii] and [N ii] lines are lower 
that that of the [Nen] 12.81/jm line, so these lines are dominated 
by emission from larger radii, where the rotation speed of the disc 
is smaller. Consequently, as in the case of a disc with a large inner 
hole, the lines do not appear double-peaked. 

So far I have neglected the effects of X-ray ionization in pro- 
ducing [Ne n] e mission, but when discu ssing observations this can- 
not be ignored. lOlassgold etaT feOOf) modelled the effects of X- 
ray irradiation in detail (see also Gorti & Hollenbach 2008). Unlike 
in the UV case, the X-ray luminosities of T Tauri stars are well- 
known, and the predicted [Ne ii] emission from X-ray ionized gas 
is comparable in strength to the emission from the photoevapora- 
tive wind (~ IO^^Lq). X-ray heating likely creates a hot, but static, 
disc atmosphere, and in real systems the emission from this atmo- 
sphere will be added to the line profiles modelled here. The line 
profiles from such a bound atmosphere have not been modelled 
in deta il, but their general shape was discussed bv lGlassgold et al] 
( l200f ). They find that emission is dominated by gas at small radii 
(< 15AU), with most of the emission arising at radii where the ro- 
tational speed of the disc gas is ~ lOkms"'. The gas temperature 
in this region is 1000-5000K, resulting in an intrinsic (Doppler) 
line- width of ~ Ikms"'. The line should therefore be double- 
peaked when viewed edge-on, with a profile similar to that pre- 
dicted here. In the face-on case, however, the profile should be 
different: narrow (~ Ikms"'), and centred on zero velocity. The 
broader (~ lOkms"') line- width predicted here is unlikely to be 
matched by a static disc atmosphere unless any turbulence in the 
disc is highly supers onic (which is generally not the case in proto- 
planetary discs, e.g.. lBalbus & Hawlevlll99a) . or unless the emis- 
sion originates very close to the star (< lAU). In addition, a static 
atmosphere results in line emission centred on zero velocity, dis- 
tinct from the blue- shifted profile of the photoevaporative wind. 
Critical to observing this blue-shift are the relative line fluxes from 
the static atmosphere and the wind. If these are comparable, as pre- 
dicted for Lx ~ 10^°erg s"' and <1> ~ 10'"s"', then the diagnostic 
blue- shift will likely be 2-5 km s~' . 

Recently iHerczeg et al.l JioOTh observed the [Neii] 12.81pm 
line from the nearby face-on (; ^7°) disc TW Hya, at a resolution 
of /l/A/i X 30,000. The emission was unresolved at an effective 
spatial resolution of ~ 40 AU, suggesting that it arises within this 
distance of the central star. iHerczeg et alj ( l2007h measured a line- 
width (FWHM) of 21 ± 4kms"', centred on -2 ± 3kms"'. The 
blue-shift is not statistically significant, in part because of the lim- 



itations of the wavelength calibration, but the line width is rather 
larger than expected from a static disc atmosphere alone. In addi- 
tion there was some evidence for asymmetry in the profile, with an 
enhanced flux on the blue side of the line, but not at a statistically 
significant level. The observed profile is consistent with emission 
from a photoevaporative wind (centred on ^ -5-lOkms ') com- 
bined with the emission from a heated disc atmosphere (centred on 
zero velocity). Further such observations are expected in the near 
future, and in this Letter I have shown that high-resolution obser- 
vations of the [Nen] 12.81yum line can provide a critical test of 
models of disc photoevaporation. Detection of blue-shifted [Ne ii] 
emission would provide unambiguous evidence of a photoevapo- 
rative wind, and observations of these line profiles may represent 
the most readily observable diagnostic of central star-driven disc 
photoevaporation. 
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